Resonance Interaction Induced by Metal Surfaces Catalyses Atom Pair Breakage 
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We present the theory for retarded resonance interaction between two identical atoms at arbitrary 
positions near a metal surface. The dipole-dipole resonance interaction force that binds isotropically 
excited atom pairs together in free space may turn repulsive close to an ideal (totally reflecting) metal 
surface. On the other hand, close to an infinitely permeable surface it may turn more attractive. We 
illustrate numerically how the dipole-dipole resonance interaction between two oxygen atoms near 
a metal surface may provide a repulsive energy of the same order of magnitude as the ground-state 
binding energy of an oxygen molecule. As a complement we also present results from density- 
functional theory. 
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Two identical, isotropically excited, atoms can be 
bound together in free space due to attractive retarded 
resonance interaction. Interestingly, Jones et al. demon- 
strated the influence of retarded resonance interaction on 
the binding^nergy of Na2 molecules [l|, [2| • Ninham and 
co-workers [3| demonstrated that, due to too drastic ap- 
proximations, the underlying theory of resonance inter- 
actions in free space derived from perturbative quantum 
electrodynamics is correct only in the non- retarded limit. 
Mechanisms of breakage and formation of molecules near 
surfaces are the key problems of catalysis. We present 
here a mechanism not previously noticed. It derives from 
the retarded resonance interaction between two identical 
atoms near a surface. Hopmeier et al. [J| provided exper- 
imental evidence for enhancement of dipole-dipole inter- 
action in microcavities and Agarwal and Gupta [5] veri- 
fied this theoretically. A way has recently been predicted 
by which very large molecules may form by jresonance 
interaction between atoms in a narrow cavity [6] . 

The aim of this brief report is to present the theory 
for retarded resonance interaction between two identical 
atoms in an excited state with the line joining the two 
atoms being perpendicular to the surface. A schematic 
illustration of the system is shown in Fig.[T] This choice 
of configuration was used because it enables us to derive 
analytical results for a configuration that has not been 
previously studied. Previous work considered only the 
case with two atoms at the same distance from the sur- 
faces. This is now extended to enable consideration of 
two atoms at arbitrary positions near a metal surface. 
We find that attractive resonance interaction between 
isotropically excited atom pairs in free space may turn re- 
pulsive close to an ideal metal surface and more attractive 
close to an infinitely permeable surface. Close to an ideal 
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metal there is a repulsive contribution which decreases 
as l/zl (where z^ is the distance between one atom 
and the surface image of the second atom) . The calcula- 
tions of resonance interaction between excited-state atom 
pairs and Casimir-Polder interaction between ground- 
state atom pairs require the proper Green's functions. 
The non-zero matrix elements of the field susceptibility 
matrix {Ta) needed for our calculations were adapted 
from the Green's functions given in the literature. Buh- 
mann and co-workers p, Q studied Casimir-Polder inter- 
action between two ground-state atoms placed with dif- 
ferent orientations near a surface. As an illustrative ex- 
ample we present numerical calculations on the resonance 
interaction between two oxygen atoms near a metal sur- 
face. The bond enthalpy for an O2 molecule is around 498 
kJ mol^^ — 5.2 eV. The sun-light spectrum has an energy 
range of 1-3 eV, so apart from solar energy one needs an 
extra 3 eV to split an oxygen molecule. A contribution to 
such bond splitting may occur due to excitation-induced 
resonance interaction between oxygen atoms near a metal 
surface. We also present results for two Zn atoms where 
the binding energy is much smaller, ^ 0.4 eV, and the 
resonance interaction effects are much stronger. 

We first recall the standard argument: Consider two 
identical atoms where one initially is in its ground state 
and the other is in an excited state. This state can also be 
represented by a superposition of states, one symmetric 
and one antisymmetric with respect to interchange of the 
atoms. While the symmetric state is likely to decay into 
two ground-state atoms, the antisymmetric state can be 
quite long-lived. The system can thus be trapped in the 
antisymmetric state [3v9]. The energy migrates back and 
forth between the two atoms until either the two atoms 
move apart or a photon is emitted away from the system. 
First order dispersion interactions are caused by this cou- 
pling of the system, i.e. the energy difference between 
the two states is separation (p) dependent. After writing 
down the equations of motion for the excited system it 
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Figure 1. (Color online) Two atoms in vacuum near a metal 
surface. Also shown are their mirror images 



is straightforward to derive the zero temperature Green's 
function for two identical and isotropic atoms |3, 10]. The 
resonance frequencies of the system are given by the fol- 
lowing equation [3| : 



I ~ a{l\uj)a{2\uj)T{p\uj)^ = 0, 



(1) 



where a{j\uj) is the polarizability of atom j [Tl|. We 
present in Fig. [5] the numerically evaluated polarizability 
of oxygen and zinc atoms used in this study. The atomic 
polarizability data were obtained from quantum chem- 
ical response theory calculations based on the complex 
polarization propagator (CPP) approach] 12] with imple- 
mentations made in the Dalton program [13] as described 
in Refs. (l4| - [la ]. In the CPP framework, damping terms 
are introduced in the response functions so as to describe 
relaxation in the system. These damping terms may be 
viewed as causing the frequency argument of the linear 
response function to become complex, and in the present 
work we employ this feature to determine the polariz- 
ability for the special case of imaginary frequencies. We 
have earlier described and applied this technique in sev- 
eral papers, see, e.g., Ref. [l7| . 

The electronic configuration of the Zn atom is 
lAr]3(i^*'4s^ and we chose to describe the electronic 
structure of this singlet ground state at the Kohn- 
Sham density functional theory (DFT) level with use 
of the Coulomb attenuated hybrid exchange-correlation 
functional CAM-B3LYP [11 [Hi in conjunction with 
Dunnin g's triple-augmented triple-^ basis set (t-aug-cc- 
pVTZ) [201 . For the O atom we consider the reference 
state of triplet spin symmetry and we describe the wave 
function by means of the multiconfigurational self- con- 
sistent field (MCSCF) method. We apply a valence active 
space, which means there is only one triplet spin-adapted 
configuration with Ms = 0. Also in this case, we adopt 
the t-aug-cc-pVTZ basis set [21 j . 

In the case of two identical atoms the above resonance 
condition can be separated in one antisymmetric and one 
symmetric part. Since the excited symmetric state has 
a much shorter life time than the antisymmetric state, 
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Figure 2. (Color online) The polarizabilities of oxygen and 
zinc atoms. The discrete Matsubara frequencies are indicated 
by the circles. Note that we have placed the n — value at 
the leil, vertical a:xes. 



the system can be trapped in an excited antisymmetric 
state [3j . The resonance interaction energy of this anti- 
symmetric state is. 



U{p) = h[uJr{p) — CJr(oo)]. 



(2) 



Since the relevant solution of Eq. ([T]) really is the pole of 
the antisymmetric part of the underlying Green's func- 
tion we can in a standard way [22] deform a contour of 
integration around this pole to obtain a both simple and 
exact expression for the resonance interaction energy, 

U{p) = _E I Jo" dujln [1 + a(l\iu;)Tu{p\iuj)] 

^""e ^ Jo°" du;a{l\^u;)Tu{p\^u;). ^^^ 

i—x,y,z 

In order to compare with previous calculations on 
dipole-dipole resonance interaction we use the approxi- 
mate linearized expression. In this way we ignore correc- 
tions due to an attractive van der Waals (Casimir-Polder) 
term. We also ignore finite size effects which will change 
the interaction for separations smaller than two atomic 
radii. 

To account for the temperature (T) dependence we 
simply replace the integration over imaginary frequencies 
by a summation over discrete Matsubara frequencies ^22| , 
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2TrkBTn/h, (4) 



where ks is the Boltzmann constant and the prime indi- 
cates that the n = term should be divided by 2. 

We consider the case when the distance, x, between 
the two atoms in the plane of the surface is zero and the 
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Figure 3. (Color online) The resonance interaction energy 
between two oxygen atoms situated one outside the other near 
an ideal metal surface. The closest atom is at 1 A from the 
surface. For comparison we have added a curve (DFT-curve) 
showing the potential for two oxygen atoms at an Au (111) 
surface as obtained from a DFT calculation including a van 
der Waals functional. Se the text for details 



distances from atoms a and b to the surface is Za and Zb, 
respectively. In other words the atoms are along a line 
perpendicular to the surface; atom a is furthest away 
from the surface (see Fig.[T]). The field susceptibility is 
a sum of three terms. Apart from the free space field 
susceptibility T^", there are p and s field susceptibility 
corrections due to the presence of a surface, T[^ and T/j, 
respectively. In free space the field susceptibility matrix 
T(p|zw) has the following non-zero matrix elements: 



Figure 4. (Color online) The resonance interaction energy 
between two oxygen atoms situated one outside the other near 
an ideal metal surface. The closest atom is 10 A from the 
surface. 



respectively. Here Ts and Tp are the reflection coefficients 
for s- and p-polarized waves, respectively, impinging on 
the surface. For an ideal metal rs — —1 and Vp — 1. For 
an infinitely permeable surface one gets a sign reversal of 
both terms. Therefore the corrections to the resonance 
interaction due to a surface have opposite sign at an ideal 
metal surface as compared to an infinitely permeable sur- 
face. The result in free space in the retarded limit is 



U{z-) « -2hca{0)/{7Tzt). 



(8) 



Near a surface we find a surface-induced nonretarded 
resonance interaction between the two atoms, 
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IUJ„ 



TU^-n) ^ -{^ + ^ + jrY- 



T^Mujn) = 2(^ + J^) --^"'-^\ 

(5) 
Here we define z^ ~ Za ^ Zb and z- = Za — Zb, with the 
first being the distance between atom a and the image 
of atom b inside the surface and the second being the 
distance between the two atoms (see Fig.[T]). The cor- 
responding surface-induced corrections to the field sus- 
ceptibility matrix [3, 18| have contributions from p and s 
polarizations (with 77 = uinZ^/c), 



U{z+) « 4:hrp/{nzl) I duja{iuj) « 2rptiuj.ja{0)/{zl,), 



(9) 

where the London approximation, a{j\iu:) ~ 
a(j|0)/ (1 -I- w^/cj^), has been used to find the fi- 
nal analytical expression. It is also possible to derive 
analytical results for the surface corrections to the 
resonance interaction at finite temperature in the 
retarded limit when a(ia;„) « a(0). Defining a;„ = nwi, 
771 = wiz_|_/c, and 7 = (—1 + e^'^)~^, we find for the p- 
and s-polarized contributions. 
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(7) 



UP{z+) « "•^'^^r^"^ {4 + 2r7h' 



+ (4 + 4771 + T]l)j + (4?7i + 3772)72} 



and 
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Figure 5. (Color online) The resonance energy between two 
oxygen (zinc) atoms both adsorbed on an ideal metal surface 
(za = zt, =1 A), solid (dashed) curve. The field susceptibility 
for this case is not described in the text. The different con- 
tributions can be obtained using the Green's functions given 
in the literature [3, Q . 



respectively. The retarded asymptote at zero tempera- 
ture then is 

U\z+) + C/''(z+) « ;ica(0)(10rp - 2r,)/(7rz^). (12) 

It is interesting to see that this asymptote has the 
same power law as the Casimir-Polder energy between a 
ground-state atom and an ideal metal surface ^23j , which 
goes as C/^^(z) w -3Sca(0)/(87rz4). 

At finite temperature the long range resonance inter- 
action driven by entropy is given by 



t/^(z+) + C/^(z+) 



8rpfci3ra(0)/z^. 



(13) 



Now let us consider the total resonance interaction in 
the retarded limit but at sufficiently small separations for 
the zero-temperature results to apply. It is given by 



U « [rica(0)A]{[(10r^ - 2r,)/(z4 )] _ 2/zl}. (14) 

It is then clear that in the case of two atoms near an 
ideal metal surface the total resonance interaction can be 
either attractive or repulsive depending on z+ and z_ . 

We will now proceed to calculate numerically the re- 
tarded resonance interaction at 300 K between two oxy- 
gen atoms near an ideal metal surface. We present in 
Figs.[3]and|4]the retarded resonance interaction between 
two atoms near an ideal metal surface. The total poten- 
tial is the sum of the contributions from free space f/o, the 
correction due to the surface from s waves [/*. and the 



correction due to the surface from -p waves [7^. In the first 
case, when one atom is adsorbed at the surface, the in- 
teraction is dominated by p-polarized surface corrections. 
When the atom closest to the surface is some distance 
away from the surface the interaction for small atom- 
atom separations is dominated by the free space contri- 
bution. However, as the distance between the atoms in- 
creases repulsive surface corrections become increasingly 
important, as can be seen in Fig. 01 In Fig. 131 we also 
compare the total potential with the potential energy for 
two oxygen atoms on a gold surface calculated within the 
Kohn-Sham DFT model with a nonlocal van der Waals 
correlation functional (rPW86/vdw-DF2) [2i,[2i|. In the 
atomistic calculation the An (111) surface is modeled by 
a 6x6x6 supercell slab structure with 44 A vacuum layer. 
The fixed oxygen atom is relaxed at the fee hollow surface 
site, and the position of the second atom is varied in the 
z direction. We find that the non-resonant contribution 
to the interaction affects the result for very small separa- 
tions of the order of 3 A or less. We have in recent work 
found that this is also the limit where finite-size correc- 
tions influence the interaction. Our results are therefore 
only meaningful when the two atoms are at separations 
beyond close contact. 

Extending the theory presented above, we furthermore 
calculated the retarded resonance interaction between 
two atoms both adsorbed on an ideal metal surface. We 
show the result in Fig. [5] as a function of the lateral dis- 
tance X between the two atoms. Surface corrections to 
the resonance interaction between atoms near an ideal 
metal surface may result in a significant repulsion. The 
DFT calculations of zinc and oxygen on the gold surface 
reveal that these atom pairs have a very weak chemical 
bond when both atoms are adsorbed on the true gold 
surface. 

We have discussed how the presence of a surface in- 
fluences the retarded dipole-dipole resonance interaction 
between two identical atoms in an isotropically excited 
state. For very small atom-atom separations, close to 
direct contact, there will be important flnite-size correc- 
tions and influence from non-resonant contributions such 
as non-retarded van der Waals forces. Beyond these lim- 
iting separations the excitation-induced interaction be- 
tween two atoms near an ideal metal can produce a re- 
pulsive energy of the same order of magnitude as the 
bonding energy of diatomic molecules. Catalytic effects 
due to the presence of a surface may induce large enough 
energies to break up or form molecules. The O2 reduction 
reaction is one of the most studied, not only because it is 
interesting scientiflcally but also important for practical 
purposes. Such catalytic effects occur, for example, at 
the cathode electrode of fuel cells, where Ft, Au, and re- 
lated alloys have been used or suggested as the catalyst. 
Yotsuhashi et al. [26] studied dissociative adsorption of 
O2 molecules on Pt and Au surfaces. They found that 
the chemical reactivity is substantially larger on Pt as 
compared to on Au surfaces. Another interesting sur- 
face is Cu(llO) where the energy gain by dissociation of 



an O2 molecule on the surface is approximately 1.3 eV 
per molecule ^STj . There are several diatomic molecules 
and metallic surfaces that will show varying strengths of 
the chemical interactions. For instance, we observe that 
the O2 molecule has much larger bond enthalpy in vac- 
uum (^5 eV) compared to the bond enthalpy of the Zn2 
molecule (^0.3 eV). Also, we regard Au to be somewhat 



more comparable to an ideal metal than for example, Pt, 
and we therefore complement the calculation of resonance 
interactions with atomistic modeling of O2 on Au(lll) 
surfaces. 
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